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Abstract 

We review the application of the spectral zeta-function to the 1- 
loop properties of quantum field theories on manifolds with boundary, 
with emphasis on Euclidean quantum gravity and quantum cosmol- 
ogy. As was shown in the literature some time ago, the only boundary 
conditions that are completely invariant under infinitesimal diffeomor- 
phisms on metric perturbations suffer from a drawback, i.e. lack of 
strong ellipticity of the resulting boundary-value problem. Neverthe- 
less, at least on the Euclidean 4-ball background, it remains possible to 
evaluate the ^(0) value, which describes in this case a universe which, 
in the limit of small 3-geometry, has vanishing probability of approach- 
ing the cosmological singularity. An assessment of this result is here 
performed, discussing its physical and mathematical implications. 
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1 Introduction 



In the Euclidean functional-integral approach to quantum gravity, one deals 
with amplitudes written formally as functional integrals over all Rieman- 
nian 4-geometries matching the boundary data on (compact) Riemannian 

3-geometries hi^ and ^^2, To take into account the gauge free- 

dom of the theory, the functional-integral measure also includes suitable ghost 
fields, described geometrically by a 1-form, hereafter denoted by = ip^dx^, 

subject to boundary conditions at {lli^hi^ and {ll2-,h2^. Although a rig- 
orous definition of the Feynman sum over all Riemannian 4-geometries with 
their topologies does not yet exist, the choice of boundary conditions still 
plays a key role to obtain an elliptic boundary-value problem, which may be 
applied to the semiclassical analysis of the quantum theory. 

In quantum cosmology, it was proposed in [2] and [3j that no boundary 

conditions should be imposed at the 3-geometry ^Si,/iij, since this might 
shrink to a point in the case of the quantum state of the universe. One would 
then have to impose suitable boundary conditions only at ^S2,/i2j, by de- 
scribing the quantum state of the universe in terms of an Euclidean functional 
integral over all compact Riemannian 4-geometries matching the boundary 

data at ^E2,/i2j- Although this approach to quantum cosmology still in- 
volves a number of formal definitions, the semiclassical evaluation of the 
corresponding wave function may be put on solid ground. The 1-loop analy- 
sis is related to mathematical and physical subjects such as cobordism theory 
(i.e., under which conditions a compact manifold is the boundary of another 
compact manifold), the geometry of compact Riemannian 4- manifolds, the 
asymptotic heat kernel, the 1-loop effective action, and the use of mixed 
boundary conditions in quantum field theory (see below). 

In particular, over the last decades many efforts have been produced to 
evaluate 1-loop quantum amplitudes for gauge fields and the gravitational 
field in the presence of boundaries, either by using the space-time covari- 
ant Schwinger-DeWitt method or the mode-by-mode analysis which relies 
on zeta-function regularization. The main motivations were the need to un- 
derstand the relation between different approaches to quantum field theories 
in the presence of boundaries and the quantization of closed cosmologies. In- 
deed, boundaries play an important role in the Feynman approach to quan- 
tum gravity as we just said, in choosing Becchi-Rouet-Stora-Tyutin- (BRST-) 
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covariant and gauge- invariant boundary conditions for quantum cosmology 
and in studying different quantization and regularization techniques in field 
theory. In particular, for the latter problem, discrepancies were found in the 
semiclassical evaluation of quantum amplitudes by using space-time covariant 
methods, where the scaling factor of 1-loop quantum amplitudes coincides 
with the Schwinger-DeWitt A2 coefficient in the heat-kernel expansion, or 
instead a mode-by-mode analysis, for which the resulting equations obeyed 
by the eigenvalues are studied through zeta-function methods [H [5]. 

If one reduces a field theory with first-class constraints [6] to its phys- 
ical degrees of freedom before quantization, one of the main problems is 
whether the resulting quantum theory is equivalent to the theories relying 
on the Faddeev-Popov gauge-averaging method or on the extended-phase- 
space Hamiltonian functional integral of Batalin, Fradkin, and Vilkovisky, 
where one takes into account ghost and gauge modes. We will see that, 
in a mode-by-mode evaluation of the covariant functional integral including 
gauge-averaging and ghost terms, after doing a 3+1 split and a Hodge-like 
decomposition of the components of metric and ghost perturbations, there 
are no exact cancellations between contributions of gauge and ghost modes, 
when linear covariant gauges are used. This lack of cancellation turns out to 
be essential to achieve agreement between different techniques. 

In [7], the C(0) calculation was performed for gravitons by restricting the 
functional-integral measure to transverse-traceless perturbations in the case 
of fiat Euclidean 4-space bounded by a 3-sphere. In [51 El [ID] this result 
was generalized to the part of the Riemannian de Sitter 4-sphere bounded 
by a 3-sphere. Both results did not coincide with those obtained by a space- 
time covariant method [11]. Hence the natural hypothesis arises that the 
possible non-cancellation of the contributions of gauge and ghost modes can 
be the cause of the discrepancy. In the work presented in [T^ [T3] such a 
suggestion was checked for the electromagnetic field on different manifolds 
and in different gauges. 

In [13] the asymptotic heat kernel for second-order elliptic operators was 
obtained in the case of pure and mixed boundary conditions in real Rie- 
mannian 4-manifolds, and in [TH] this analysis was improved. In the light of 
these results, the conformal anomalies on Einstein spaces with boundaries 
were recalculated in [T6] . 

In [T7] the linearized gravitational field was studied in the geometric 
framework of [7| (i.e. flat Euclidean 4-space bounded by a 3-sphere), and 
the resulting ({0) value was compared with the space-time covariant calcula- 
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tion of the same Faddeev-Popov amplitudes, by using the corrected geomet- 
ric formulae for the asymptotic heat kernel in the case of mixed boundary 
conditions [16]. 

However, in the case of mixed boundary conditions involving tangential 
derivatives of metric perturbations, no geometric formulae for the asymptotic 
heat kernel are available as yet, and one has to resort, to the best of our 
knowledge, to analytic techniques along the lines of the work in [TSl [TU] . 
This is what our review is mainly devoted to. For this purpose, section 2 
derives the integral representation of the spectral zeta-f unctions, as obtained 
in [8l|20l[2T]. Sections 3 and 4 are devoted to C(0) values for scalar and gauge 
fields, respectively. Detailed calculations for the gravitational field begin in 
section 5 and continue until section 7. The strong ellipticity issue is studied 
in section 8, with examples. Concluding remarks and open problems are 
presented in section 9. 

2 Integral representation of the spectral zeta- 
function 

A convenient method for the calculation of the spectral zeta function for the 
case when the spectrum is not known explicitly, but only the structure of 
the basis functions of the corresponding differential operator are known, was 
proposed in |8]. Here we sketch the basic ideas and formulae of this method; 
for a related approach see Dowker |22l |23] . 

Let us consider the second-order operator F, which represents the second 
functional derivative of the Euclidean action of the model under consideration 
with respect to the field variables. It is convenient to single out the mass 
term from the operator F . As a manifold we consider the part of the 
closed Euclidean de Sitter space ("Euclidean ball"). Then suppose that we 
have a full set of basis functions u\{T\rn?) of this massive operator, i.e. 



Here r is the Euclidean time parameter, lower case Latin indices i enumerate 
the modes of the field variables while capital Latin indices enumerate the 
basis functions, the coefficients of the second-order derivatives with 

respect to the time parameter. 




(1) 
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The only condition which these basis functions should satisfy is regularity 
in the Euclidean ball < r < r_|_. Then the eigenvalues A of the operator 
F{d/dT) + m^a with homogeneous Dirichlet boundary conditions satisfy the 
equation 

u'Ar+\m'-\) = 0. (2) 

For other types of boundary conditions the basis functions u in Eq. ([2]) should 
be substituted by the corresponding combination of basis functions and their 
derivatives. For the case of Dirichlet boundary conditions, which we consider 
in more detail in this section, the equations defining all the eigenvalues can 
be collected in one equation, 

detu^(r+|m2 - A) = 0, (3) 

in which the determinant is taken with respect to the indices i and A of 
the square matrix u\. Then using the Cauchy formula and the well-known 
relation between det and tr, we can rewrite ({s) as an integral, 

over the contour C in the complex plane of z, which encircles all roots of ([3]). 

It is necessary to note that positivity and real-valuedness of roots of ([3]) 
is guaranteed by self-adjointness and positive definiteness of F^d/dr) +'w?a, 
which is assumed here. 

It should be stressed also that each basis function may be taken with an 
overall normalization factor depending on {m? — z). This can lead to the 
additional roots of Eq. ([3]), which are irrelevant to the eigenvalues of the 
elliptic operator under consideration. To avoid such an effect, consider the 
singular point of the radial equation at r = 0. The asymptotic behaviour of 
u\{T\m? — z) for r — )■ has, according to the asymptotic expansion theory, 
a power-law form [2^ 

u{T\m^ - z) ^ UqT^ + 0{t^+^), r^O, (5) 

where is a positive integer number. Thus, to avoid the additional roots of 
Eq. ([3]), it is sufficient to require that uq should be independent of the mass. 

If we now assume that the basis functions are analytic in the complex 
plane of the mass variable m^, then we can continuously deform the original 
contour of integration C to the new contour C*, which encircles the cut in the 
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complex plane of the functions , coinciding with the negative real axis. 
Thus, the general expression for C(s), to be analytically continued to s = 0, 
looks like 

For the analytic continuation of (EI) from the convergence region domain to 
s = 0, take into account that the contour C includes the two boundaries 
of the negative real axis and a circle around the point 2; = of some small 
radius e. Therefore 

, sin(7rs) dM^ d , , , o ,.2^ 

where the first term is a jump of the integrand in ([6]) on the cut of the 
function , integrated along this cut z = — M^. 

Let us transform Eq. ([7]) by the following sequence of operations: first 
analytically continue both terms into the neighborhood of s = and then go 
to the limit e = 0. The integral along will vanish because of the regularity 
of M(r+|m^ — z) at z = 0. 

It is not so difficult to show that for a quantum-mechanical system with 
a finite number of degrees of freedom, as s — )■ we have 

as) = Iio, + s[I]^ + Ois'), (8) 

where 

I{M^) =ti\nu{T+\m^ + M^), (9) 

/log is the coefficient of InM^ in the expansion of / as — )■ 00, [J]°° is the 
regular part of this expansion as — )■ 00 and [/]o = /(O), [I]"^ = [/]°° — [/Jo- 
lt is obvious that in this case 

C(o) = /log, (10) 
C'(o) = [/]^. (11) 

This result shows that C(0) and C'(0) get a contribution from the asymptotic 
value of the basis fucntion u\{t\M'^) for — )■ 00. But their asymptotic 
behaviour can be obtained from the JWKB approximation for the corre- 
sponding equation [2^ . 
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The problem becomes much more comphcated when we study field theo- 
ries, which have an infinite number of modes, because the trace in (jS]) becomes 
the divergent series 



The question arises of how the parameter s can regularize this divergent 
series. Let us interchange the summation and integration operations in (fT2|) . 

and consider the asymptotic behaviour of the integral 

f^^J\nuir^\m^-z)]i, (14) 

for the collective index A growing to infinity. The numerical parameter tend- 
ing to infinity with the growth of A is the parameter n enumerating the har- 
monics of the radial equation. Thus, the question of convergence for the sum 
f[T^ reduces to the analysis of the asymptotic behaviour of (HM for n — )■ cxd. 
Fortunately, the so-called uniform JWKB expansion for the basis functions 
has an important property [211 [251 when it is considered as a function of 
the two arguments n — )■ oo and the ratio z/n^, 

\iau{T\m^ - z) = v?jwKB (n^, ^) , (15) 

then it is uniform in the second argument, < \z\/n^ < oo, and at most 
has a power-law growth of finite order k in the first argument, — > oo. 
Therefore, substituting f[T5|) into f[T^ and making the change of integration 
variable z — )■ n'^z one finds that this integral has an asymptotic behaviour, 

which converges for some s > and, due to uniformity of v'jwkb(^^, z) in z, 
has a bound const ^{ri^)^"^ providing the convergence of the infinite series 
f[T3|) for some large positive s. Thus, large values of s > regularize the 
divergent sum in f lT^ . 
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Making the change of integration variable z — )■ n^z in Eq. f ll2p and 
interchanging back the order of integration and summation one can represent 
the ( function in the form 

where I{—z, s) is the manisfestly regularized infinite sum 

^(-^'^) = E4j[1^«(^+I^'-^)]a- (18) 

A ^ 

Similarly to (JTj), we can split the integral ( ITTll over C into a sum of 
two terms and show that the integral around the circle Ce tends to zero. 
However, unlike models with a finite number of physical variables, the series 
([T8|l analytically continued from its convergence domain generally has a pole 
at s = 0, 

I{M^,s) = - ^ + /^(M2) + 0(s). (19) 

s 

Therefore, instead of the formula (IHD we obtain the following result for the 
case of field theory: 



+ 1^ 



r-Rioo 
Jo 



fdMHnM^^^^^^yOis'), (20) 



where ()iog and []^ have the same sense as in the quantum-mechanical case. 
Thus we have 

C(o) = (/^)iog + [P°'T^ (21) 

C'(0) = [IX - lyMHnAP^^^^^^. (22) 



These equations generalize the algorithms f flQl) and ( JTTj) to field theories 
with an infinite number of physical modes. But these generalizations are 
non-trivial: only the terms (/^)iog in f l2T|) and [I^]"^ in f l22|) are similar to 
the expressions f lTOj) and ffTTj) . The terms including 1^°^'^ do not have analogs 
in a theory with a finite number of modes. These terms are responsible for 
the non-trivial renormalization of the ultraviolet divergences performed by 
the ^-function regularization. 
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The formalism just described is a fine tuned scheme for the computation 
of C(0) and C'(0)- It does, however, not allow to extract other properties 
of C{^)- But ill order to determine heat kernel coefficients or the Casimir 
energy associated with some quantum field theory models, other particular 
properties are needed [26l |271 [28] . These can be found by analytically contin- 
uing the zeta function ({s) as given in Eq. f[T^ to a meromorphic function 
in the complex plane. The details of this procedure depend very much on 
the explicit form of M(r+|m^ — z) in Eq. (fT3|) . In general one can only say 
that adding and subtracting the asymptotic expansion briefiy outlined in 
Eq. (fT5l) is crucial to the method, but the precise nature of integrals and 
series to be done to obtain the analytical continuation depend on exactly 
what ifjwKBin'^iz/'n?) actually is. For the example of the scalar Laplacian 
on a four dimensional ball with various boundary conditions details will be 
provided in the next section. 

3 Dirichlet, Neumann and Robin Boundary 
Conditions 

One becomes familiar with Dirichlet boundary conditions as soon as one 
studies potential theory. The first boundary-value problem of potential the- 
ory is the existence of a function, harmonic in a closed region, and taking 
on preassigned continuous boundary values. This is known as the Dirichlet 
problem, and is the oldest existence theorem in potential theory. Usually, 
one first tries to express a harmonic function in terms of its boundary val- 
ues. One then sees if the expression found continues to represent a harmonic 
function when the boundary values are any given continuous function. 

The problem of finding a function, harmonic in a region, and having 
normal derivatives equal to the function given on the boundary is instead 
the Neumann problem, or the second boundary-value problem of potential 
theory. The theorem asserting the existence of a solution of this problem is 
known as the second fundamental existence theorem of potential theory 

In the semiclassical approximation of the quantum theory of a real scalar 
field in a real Riemannian background with boundary, the guiding principle 
for the choice of boundary conditions is that the boundary data should refiect 
those particular conditions which lead to a well-posed classical boundary- 
value problem. Thus, on using the background-field method, the scalar-field 
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perturbations ip are required to obey one of the following three boundary 
conditions on the bounding surfaces 

(i) Dirichlet problem: 

(ii) Neumann problem: 



(/? = at dM, (23) 



^ = at 9M, (24) 

OT 



(iii) Robin problem: 



^ + -<^ = at dM. (25) 
or T 

For example, in the case of a massless scalar field at 1 loop about flat 4- 
dimensional Euclidean space bounded by a 3-sphere, the technique of section 
2 may be used to find the following values for the resulting anomalous scaling 
factors: 

Cd(0) = -^ (Dirichlet), (26) 
29 

Cn{0) = — (Neumann), (27) 

Ci?(0) = -j^ - ^(n - 1)3 (Robin). (28) 

In this particular case, the ({0) values coincide with the conformal anomaly, 
since massless scalar field theories are conformally invariant in fiat space-time. 
It was not until in [21] that a powerful analytic algorithm was developed for 
the analysis of the Robin case, and the first correct geometric results for 
C(0) were only published in [32] and [T5| 116] . More recent work on real 
scalar fields on the Euclidean ball in various dimensions can be found in 
[201 EH [221 ESI [33]. 

In order to outline the contour integration method for the analysis of 
spectral zeta functions as a function of s let us exploit this opportunity 
and rederive Eq. (126|) with an indication on how to obtain Eqs. (1271) -( |28|) . 
Given the treatment of arbitrary dimension does not cause any additional 
complications, we will consider the D = d + 1 dimensional ball [2S] • 

A massless scalar field leads to the eigenvalue problem for a Laplacian, and 
for the spherically symmetric problem at hand the use of polar coordinates 
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(29) 



seems appropriate. In these coordinates the Laplacian reads 

92 d d 1 , 

with A^v" the Laplacian on the d-dimensional sphere, = S*^. 

By imposing Dirichlet boundary conditions on the sphere, the boundary 
of the ball, eigenvalues are determined by the transcendental equation 

Ju{\) = 0, (30) 

with u = £ + {d — l)/2, £ = 0, 1,2, and with the radius R of the ball being 
chosen as i? = 1. The degeneracy d^ of each eigenvalue equals the degeneracy 
of the eigenvalues of the Laplacian on the cZ-sphere and for d > 2, the case 
we will concentrate on in the following, it equals 

{i + d-2)\ 



d^ = (2£ + d- 1)- 



i\{d-l)\ 



(31) 



The information provided suffices to give an explicit representation of the 
associated zeta function as given in Eq. ( |T3|) . i.e. 



((■■>) 



c 



sm TVS 

71 



(32) 



where this equality is obtained by deforming the contour C to the imaginary 
axis. The relevant uniform asymptotic behaviour outlined in Eq. (fT5|) for 
the given example follows from 



V2^ii + z'y/' 

k/u is fixed [M 13 



Uk{t) 



k=l 

Here t 



(33) 



and 



valid for u ^ 00 as z 
1] = yl + z^ _|_ ln[z/(l + a/1 + 2;2)]. Higher powers in t follow from the 
recursion 



Uk+l(t) 



t 

\t\l-e)u',{t)+^- I rfr(l-5r>,(r). 
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starting with 'Uo('^) = 1- On defining polynomials Dn{t) from the expansion 



In 



k=l 



n=l 



Dnit) 



(34) 



the leading few polynomials are 



D2{t) 



-t f, 

8 24 ' 

1 2 3 , 

—t^ - -e 

16 8 



25 



531 



16 



(35) 



221 



1105 



384 640 128 1152 
with many more polynomials easily found by using an algebraic computer 
program. 

By adding and subtracting leading terms, for (^(0) in D = 4 we will 
ultimately choose = 3, the zeta function (132|) splits into the pieces (after 
substituting k = zu) 

N 

Cis) = Z{s) + J2 Ms), 



i=-l 



where 



Z{s) 



sin(7rs) 



TT 



^2 



-2s 



-In 



V2^ii + z^y/^ 



d_ 

dz 

N 



I In " Iu{zu 



(36) 



n=l 



and the Ai{s) result from the different orders in the asymptotic expansion, 
explicitly 



Aois) 
AAs) 



sin(7rs) 



TT 



Sm TTS 



TT 



Sm TTS 



TT 



oo 

^ c/, y dz{ziy)-^'-^\n{z-''e'''^) , 



di, / dz{ 



zu 



'2s 



d_ 

dz 



Hi + z')-'i\ 



d fD.,{t) 



dz \ V 
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It can be shown that Z{s) is analytic in the half-plane {d — 1 — N)/2 < 3f?s. 
These formulas therefore make it possible to find a representation of ({s) 
valid for any value of s. 

Choosing suitably large, given the factor sin(7rs) in f l36|l . Z{s) will 
therefore not contribute to C(0)- As far as ({0) is concerned, it therefore 
suffices to only consider the Ai{s) further. By introducing the so-called base 
zeta function, 

Cu{s) = J2duiy-'% (37) 
y4_i(s) and Aq{s) are readily evaluated as 

A^{s) = --Cm{s). 

In order to compute the higher Ai{s) note that the polynomials -Dj(t) can be 
written as 



i+2b 

6=0 



with the coefficients Xi^b easily determined from the definition (l34l) of Di(t)] 
to read off the numbers Xi^b for i = 1, 2, 3, see also The 2;-integrals are 

then easily done and one finds 



r(s + 6 + i) 

6=0 

Concentrating on ^(0) in four dimensions we note that 



£=0 
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which allows us to write 



4A/7r r(s + 1) 



At s = we compute 



^-i(0) = (--] Cr(-3) = — 

^ ' 40F V 2; ' 240' 

Ao(0) = Ai(0) = ^2(0) =0, 

^3(0) = -^E^3,6 = -3^ 



and thus as stated 

C(o) 



2^ 720' 

6=0 



1 



240 720 180 



In the same manner Neumann and Robin boundary conditions can be treated 
by starting with the implicit eigenvalue equation 

+ J,(A) + AJ:(A) = 0, 



where u = corresponds to Neumann boundary conditions; for details about 
the very similar calculations we refer to [26j . 

For complex scalar fields, the most general case corresponds to mixed 
boundary conditions, i.e. when the real part obeys Dirichlet conditions and 
the imaginary part obeys Neumann conditions, or the other way around. In 
the light of fl26|) - fl27|) . the resulting conformal anomaly for a complex massless 
field on the Euclidean ball is found to be 

C(0) = ^. (38) 
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4 Mixed boundary conditions for gauge fields 



We are interested in the 1-loop amplitudes of vacuum Maxwell theory in 
the presence of boundaries. Since in the classical theory the potential is 
subject to the gauge transformations 



(39) 



this gauge freedom is reflected in the quantum theory by a ghost 0-form, i.e. 
an anticommuting, complex scalar fleld, hereafter denoted again by (p. The 
two sets of mixed boundary conditions consistent with gauge invariance and 
Becchi-Rouet-Stora-Tyutin (hereafter BRST) symmetry are magnetic, i.e. 



A, 



0, 



dM 



J dM 



dM 



or electric, i.e. 



^0 

dAk 
dr 



dM 



dM 



(40) 

(41) 
(42) 

(43) 
(44) 

(45) 



J dM 



where $ is an arbitrary gauge-averaging functional deflned on the space of 
connection 1-forms A^dx'^. Note that the boundary condition (H2!) ensures 
the gauge invariance of the boundary conditions fH0l) - (l4T!) on making the 
gauge transformation (139|) . Similarly, the boundary condition fH5|) ensures 
the gauge invariance of on transforming the potential as in fl39p . 

For example, when the Lorenz gauge-averaging functional is chosen. 



the boundary condition 
tensor of the boundary) 



T]) reduces to {K being the extrinsic-curvature 
dAo 



dr 



+ AqTt K 



0, 



dM 
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by virtue of Eq. f l4Up that sets to zero at the boundary all longitudinal and 
transverse modes. 

It is also instructive to prove the BRST invariance of our boundary condi- 
tions. For this purpose, take e.g. the electric boundary conditions fH3|) - fH5|) . 
jointly with the BRST transformations (the ghost 0-form corresponding to 
independent and real-valued ghost fields denoted by uj and ip, while 5\ is an 
anti-commuting gauge parameter) 

^BRST^M = (V^^)5A, (46) 

5brst^^ = (VM^)5A, (47) 
^BRST^ = 0. (48) 
Now if ip obeys Neumann boundary conditions, 



0, (49) 



J dM 

then by virtue of P6|) one finds (hereafter 6 = 5brst) 

5 (n^A'^) = n, (m^) = {5X)n, (V^V-) , (50) 

and this variation vanishes at the boundary by virtue of f l49p . Thus, the 
boundary condition 

n^^^H = 0, (51) 

which is the covariant form of ( H3|) . is preserved under the action of BRST 
transformations. Further details can be found in [IT]. 

For a given choice of one of these two sets of mixed boundary condi- 
tions, different choices of background 4-geometry, boundary 3-geometry and 
gauge-averaging functional lead to a number of interesting results. We here 
summarize them in the case of a background given by flat Euclidean 4-space 
bounded by one 3-sphere (i.e. the disk) or by two concentric 3-spheres (i.e. 
the ring). 

(i) The operator matrix acting on normal and longitudinal modes of the 
potential can be put in diagonal form for all relativistic gauge conditions 
which can be expressed as 

$,(A) = V^A^-6AoTrir, (52) 
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where denotes covariant differentiation with respect to the Levi-Civita 
connection of the background, and 6 is a dimensionless parameter. 

(ii) In the case of the disk, the Lorenz gauge (set 6 = in ( 152|) ) leads to a 
^(0) value 

SI 

a(0) = --, (53) 

for both magnetic and electric boundary conditions, which agrees [16] with 
the geometric theory of the asymptotic heat kernel. However, the C(0) value 
depends on the gauge condition, and unless b vanishes it also depends on the 
boundary conditions. 

(iii) In the case of the ring, one finds 

C(0) = 0, (54) 

for all gauge conditions, independently of boundary conditions. This result 
agrees with the geometric formulae for the heat kernel, since volume (i.e. 
interior) contributions to ({0) vanish in a flat background, while surface (i.e. 
boundary) contributions cancel each other. 

(iv) In the case of boundary 3-geometries given by one or two 3-spheres, the 
most general gauge-averaging functional takes the form [131 ES] 

^A) = 71 (^V% + Tr K - 73 ^'W'A, (55) 

where 71, 72 and 73 are arbitrary dimensionless parameters, which give differ- 
ent "weight" to the various terms in the 3+1 decomposition of V^A^. Thus, 
unless 7i,72 and 73 take some special values (cf. ( 152|) ). it is not possible to 
diagonalize the operator matrix acting on normal and longitudinal modes of 
the potential. 

(v) The contributions to C(0) resulting from normal and longitudinal modes 
do not cancel in general the contribution of ghost modes. Thus, transverse 
modes do not provide the only surviving contribution to 1-loop amplitudes. 
In other words, all perturbative modes are necessary to recover the correct 
form of 1-loop semiclassical amplitudes. 
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5 Boundary conditions for the gravitational 
field 



For gauge fields and gravitation, the boundary conditions are mixed in that 
some components of the field (more precisely, a 1-form or a symmetric tensor 
of type (0, 2)) obey a set of boundary conditions, and the remaining part of 
the field obeys another set of boundary conditions. Moreover, the bound- 
ary conditions are invariant under local gauge transformations provided that 
suitable boundary conditions are imposed on the corresponding ghost 0-form 
or 1-form. 

We are here interested in the derivation of mixed boundary conditions for 
Euclidean quantum gravity. The knowledge of the classical variational prob- 
lem, and the principle of gauge invariance, are enough to lead to a highly 
non-trivial quantum boundary-value problem. Indeed, it is by now well- 
known that, if one fixes the 3-metric at the boundary in general relativity, 
the corresponding variational problem is well-posed and leads to the Ein- 
stein equations, providing the Einstein-Hilbert action is supplemented by a 
boundary term whose integrand is proportional to the trace of the second 
fundamental form. In the corresponding quantum boundary-value problem, 
which is relevant for the 1-loop approximation in quantum gravity, the per- 
turbations hij of the induced 3-metric are set to zero at the boundary. More- 
over, the whole set of metric perturbations /i^y are subject to the so-called 
infinitesimal gauge transformations 



where V is the Levi-Civita connection of the background 4-geometry with 
metric g, and tp^dx^ is the ghost 1-form. In geometric language, the difference 
between h^i, and h^^ is given by the Lie derivative along ip of the 4- metric g. 
For problems with boundary, Eq. implies that 



where the stroke denotes, as usual, 3-dimensional covariant differentiation 
tangentially with respect to the intrinsic Levi-Civita connection of the bound- 
ary, while Kij is the extrinsic-curvature tensor of the boundary. Of course, 
ifQ and (p>i are the normal and tangential components of the ghost 1-form, 
respectively. Note that boundaries make it necessary to perform a 3-1-1 split 



(56) 



hij = hij + + Kijipo, 



(57) 
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of space-time geometry and physical fields. As such, they introduce non- 
covariant elements in the analysis of problems relevant for quantum gravity. 
This seems to be an unavoidable feature, although the boundary conditions 
may be written in a covariant way. 

In light of ( 157|) . the boundary conditions 



hi 



are gauge invariant, i.e. 



h 



dM 



dM 







0, 



(58) 



(59) 



if and only if the whole ghost 1-form obeys homogeneous Dirichlet conditions, 
so that 



dM 



dM 



0, 



0. 



(60) 
(61) 



The conditions (!60!) and (!6T|) are necessary and sufficient since ipo and ipi are 
independent, and 3-dimensional covariant differentiation commutes with the 
operation of restriction at the boundary. Indeed, we are assuming that the 
boundary is smooth and not totally geodesic, i.e. Kij ^ 0. However, at those 
points of dM where the extrinsic- curvature tensor vanishes, the condition 
fl60|) is no longer necessary |35] . 

The problem now arises to impose boundary conditions on the remain- 
ing set of metric perturbations. The key point is to make sure that the 
invariance of such boundary conditions under the infinitesimal transforma- 
tions (jSnD is again guaranteed by fpl - flMl) . since otherwise one would obtain 
incompatible sets of boundary conditions on the ghost 1-form. Indeed, on 
using the Faddeev-Popov formalism for the amplitudes of quantum gravity, 
it is necessary to use a gauge- averaging term in the Euclidean action, of the 
form 

^g.a. = t^^tV / ^'g.^^'^V^g d'x, (62) 
oZnUa J j^j 

where $^ is any relativistic gauge-averaging functional which leads to self- 
adjoint elliptic operators on metric and ghost perturbations. One then finds 
that (here <l>^ = g^,^^") 



S^^ih) = <I>^(/i) - <I>^(/i) 



(63) 
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where is an elliptic operator that acts linearly on the ghost 1-form. Thus, 
if one imposes the boundary conditions 



dM 



0, 



0, 



J dM 



(64) 
(65) 



their invariance under fl56p is guaranteed when fl6UI) and 0611) hold, by virtue 
of (1^. Hence one also has 



dM 



0, 



0. 



J dM 



(66) 
(67) 



In section 7 we shall study this scheme, first proposed in [SB], when the 
linear covariant gauge of the de Bonder type is chosen. We will see that this 
leads to boundary conditions which involve normal and tangential deriva- 
tives of normal components of metric perturbations, and the resulting 1-loop 
divergence will be evaluated. 



6 Equations for basis functions and their so- 
lutions for pure gravity 

For the reasons described in the introduction, we study pure gravity at 1 loop 
about a flat Euclidean background with two concentric 3-sphere boundaries, 
and eventually let one of the 3-spheres shrink to a point. Our approach to 
quantization follows the Feynman-DeWitt-Faddeev-Popov formalism [3T 
Hence we deal with quantum amplitudes of the form 



Z [boundary data] 



/ilM /i2M exp(-j£;) 



c 



With our notation, C is the set of all Riemannian 4-geometries matching 
the boundary data, /xi is a suitable measure on the space of metrics, /i2 is a 
suitable measure for ghosts, is an arbitrary gauge-averaging functional, 
and the total Euclidean action reads (in c = 1 units) 



/ Wi^v^d^ d'a: + -^ [ Tr iT^dit^ 
WnG Jm SttG J^m 



X 



+ 



167rG ./a/ 2a 



(68) 
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where ^^^i? is the trace of the 4-dimensional Ricci tensor. Of course, K is the 
extrinsic-curvature tensor of the boundary, q is the induced 3-metric of 9M, 
and q; is a positive dimensionless parameter. The ghost action Jgh depends 
on the specific form of Denoting by h^v the perturbation about the 
background 4-metric g^^i^, one thus finds field equations of the kind 

where is the 4-dimensional elliptic operator corresponding to the form 
of = Qufj.^^ one is working with. Here we choose the de Bonder gauge- 
averaging functional 



where is covariant differentiation with respect to and h = gf^" h^i,. 
The corresponding operator is the one obtained by analytic contin- 

uation of the standard D'Alembert operator, hereafter denoted by Q . The 
resulting eigenvalue equation is 

Now we can make the 3+1 decomposition of our background 4-geometry 
and expand /ioo, /ioi and hij in hyperspherical harmonics as 



n=l 



hoi{x,T) = ^ 



n=2 



(69) 
(70) 



n=3 
oo 



+ ^g(")(x) 



n=l 



+ 



(n^-l) ' 3 
E fn{r){sl;^{x)+S^^{x))+Ur)G^\x) 



n=3 



(71) 



Here Q W (x) , Sf' (x) and G^;^ (x) are scalar, transverse vector and transverse- 
traceless tensor hyperspherical harmonics, respectively, on a unit 3-sphere 
with metric Qj. 
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The insertion of the expansions fl69p - fl7ip into Eq. 0681) leads to the fol- 
lowing system of equations, 

l„a„(r) + Bnbnir) + C„e„(r) = 0, (72) 

5„6„(r) + ^„a„(r) + F„c/„(r) + G„e„(r) = 0, (73) 

L„rf„(r)+M„6„(r) = 0, (74) 

Nnen{r) + Pnbn{r) + Qnan{r) = 0, (75) 

HnCnir) + KJn{r) = 0, (76) 

Rnfnir) + SnCnir) = 0, (77) 

fnKir) = 0. (78) 

Since our background is flat, after setting a = 1 in fl68|) the operators ap- 
pearing in Eqs. f l72]) -f l75]) take the form (for all integer n > 3) 

S,. = 4, (80) 

Cn = ^, (81) 

^ = -(^2-1), (83) 
r 

^ _ rf2 1 (n^ + 5) , , 

^"=d^~7d7"^^ + ^"' ^^^^ 
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Qn = 6, (92) 

S„ s -, (94) 

r 

^ 1 d (n^ - 1) , 

T„ = — - -— - ^ ^ + A„. 95 

Inserting the operator T„ from Eq. (!95|) into Eq. ( I78|) one can easily find the 
basis function describing the transverse-traceless symmetric tensor harmonics 
which usually are treated as physical degrees of freedom: 

knij) = aiTln{MT) + a2TKn{MT) , n = 3, . . . , (96) 

where M = a/— A and / and K are modified Bessel functions. 

However, the equations (172l)- (!75l) for scalar-type gravitational perturba- 
tions lead to a rather complicated entangled system as well as Eqs. (I76|) and 
fl77|) . describing vector perturbations. In Refs. [121 [13], where the analogous 
problem was studied for the electromagnetic field, a method was used to de- 
couple a similar entangled system for normal and longitudinal components of 
the 4-vector potential. The idea is that one can diagonalize a 2 x 2 operator 
matrix after multiplying it by two functional matrices. In some cases one can 
choose these functional matrices in such a way that the transformed opera- 
tor matrix is diagonal and the corresponding differential equations for basis 
functions are decoupled. However, in the case of scalar-type gravitational 
perturbations we have a 4 x 4 operator matrix. To diagonalize such a matrix 
it is necessary to solve a system of 24 second-order algebraic equations with 
24 variables. This problem seems a rather cumbersome one and we thus use 
another method. For this purpose, we assume that the solution of the system 
of equations ( 172|) -( 175l) is some set of modified Bessel functions with unknown 
index u. Let us look for a solution of this system in the form 
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bn{r) = /32W,{Mt), (98) 

dn{r) = 15^tW,{Mt), (99) 

e„(r) = I3^tW,{Mt). (100) 

Here, Wy is a linear combination of modified Bessel functions 1^ and Ky 
obeying the Bessel equation 

-r^ + —, o-M^] Wy{MT) = 0. 101 



Now, inserting the functions (!97|) - fll00l) and the corresponding operators into 
the system of equations (1721) - (1751) . and taking into account the Bessel equa- 
tion (llOip . one finds the following system of equations for /32, /Ss and (^4, 



(z/2 -n^- 6)/3i + 4/32 + 2(3^ = 0, 

6{n^ - l)/3i + 3(z/2 -n^- A)(52 + 4(^2 - A)I3^ - 2(n' - l)/34 

4/32 + (z/'-n2 + 4)/33 = 0, 



6/3i - 4/32 + (z^' 



n 



2)/34 = 0. 



(102) 

(103) 
(104) 
(105) 



The condition for the existence of non-trivial solutions of the linear homoge- 
neous system (I102p - fll05l) is the vanishing of its determinant, i.e. 



f 2 2\2 



0. 



(106) 



The roots of Eq. fll06p are 



2 2 
n , V 



[n 



2)\ 



The positive values of v provide the orders of modified Bessel functions. Now 
we can write down the /3's corresponding to different values for z/'s. For v = n 
one has 

/34 = 3/3i, /32 = /33 = 0, (107) 



or 



/3i = 0, /33 = -/32,/34 = -2/32. 



For V = n — 2 one has 



(108) 



in + 1) 

/32 = (n + l)/3i, /33 = ^-^/3i, /34 = -/3i. 



(109) 
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Last, for u = n + 2 one has 



(n — 1) 



(110) 



Having the Eqs. fll07l) - flll0p we can get the basis functions for scalar-type 
gravitational perturbations (l97 l) -( fT00|) . 



an(r) = - 7i4(Mr)+73J„_2(Mr)+74/„+2(Mr) 



+ 5iKn{MT) + 5^Kn^2{MT) + 5^Kn+2{Mr) , (111) 



hn{T) = 72/„(Mr) + (n + l)73/„_2(Mr) 
- (n - l)74/„+2(Mr) + 52i^n(Mr) 
+ (n + l)53i^„_2(Mr)-(n-l)54i^„+2(Mr) 



(112) 



+ 



r 

(n-1) 



+ 11 

72/n(Mr) + ^ ^73/„_2(Mr) 



n + 2) 

(^ + 1) 
(n-2) 



n - 2) 

74^+2 (Mr) - 62K^{Mt) 



in — 1) 

53i^n-2(Mr) + ^^54ir„+2(Mr) 
(n + 2) 



(113) 



r [37i/„(Mr) - 272/„(Mr) - 73/„_2(Mr) 
74/„+2 (Mr) + 35iir„ (Mr) - 252^^™ (Mr) 
53i^„_2(Mr)-54i^„+2(Mr) 



(114) 



We can find the basis functions for vectorlike gravitational perturbations 
in a similar way. Let us suppose that 



and 



Cn{T)=e,W,{MT), 

fn{r)=e2TW,{MT). 



(115) 
(116) 
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Inserting (11151) and (11161) into Eqs. (17^ and (1771) one has the system 

-n^ - h)ei + 2{n^ - 4)^2 = 0, 

2ei + -n^ + 3)52 = 0. 
The determinant of the system (I117P is 

- n^f - 2(z/2 - n^) - W + 1, 

and its positive roots are nil. For v = n + 1 one has 

1 



£2 



(n + 2) 



^1, 



(117) 
(118) 



and for v = n — 1 one has 



£2 



in -2) 

and correspondingly the basis functions ( IllSp and ( 1116^ take the form 
c„(r) = £iJ„+i(Mr) + £2/n-i(Mr) + r/iif„+i(Mr) + r^^Kn-AMr), (119) 



/n(r) 



r 



^ £iJ„+i(Mr) + -A_£2/n-i(Mr) 



ri + 2 



(n + 2) 



7]iKn+i{Mr) 



(n-2) 
1 



fn-2: 



7]2Kn-l{MT) 



(120) 



We have also to find the basis functions for ghosts. The eigenvalue equa- 
tions for ghosts in the de Bonder gauge have the form 

and the corresponding fields can be expanded on a family of 3-spheres as 



(fi{x,T) = 

n=2 



^o(x,r) = 5;;/„(r)QW(a;), 

n=l 

n^nij)— -T^Vn[T)S\ \X) 

[n^ — 1) 
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;i2i) 



;i22) 



The functions ln{T),mn{T) and Pn{T) can be found similarly to those for 
harmonics of gravitational perturbations. They have the form 



ln{r) = - 



Kjn+l{MT) + K2ln-l{Mr)+e^Kn+l{MT)+e2Kn-l{MT) , (123) 



m„(r) = -(n- l)/ti/„+i(Mr) + (n + l)K2/n-i(Mr) 

- {n-l)e,Kn+i{MT) + {n + l)e2Kr,^,{MT), (124) 

p„(r) = ^IniMr) + pKaiMr). (125) 



7 Barvinsky boundary conditions 

As we know from section 5, one can set to zero at the boundary the gauge- 
averaging functional, the whole ghost 1-form, and the perturbation of the 
induced 3-metric. With the notation of section 5, after making an infinites- 
imal gauge transformation of the metric perturbation h^y according to the 
law (!56|) . one finds in the de Bonder gauge (cf. (!63l) ) 



$f (/.) 



□ 



(126) 



where R^^ is the Ricci tensor of the background, and the elliptic operator 
on the right-hand side of fll26p acts linearly on the ghost 1-form. In our flat 
Euclidean background, the Ricci tensor vanishes, and on making a 3-1-1 split 
of the de Bonder functional (^'fP and of the ghost 1-form (^^, the boundary 
conditions proposed in [36] read as (unlike section 5, we here consider only the 
de Bonder gauge-averaging functional, with the corresponding superscript 
dD) 



hij 



dM 



dM 



h 



dM 



0, 



J dM 



dM 



dM 



dM 



0, 



0. 



dM 



(127) 
(128) 
(129) 
(130) 
(131) 
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Once again, the vanishing of the whole ghost 1-form at the boundary ensures 
the invariance of the boundary conditions (11271) under the transformations 
fl56|l . At that stage, the only remaining set of boundary conditions on met- 
ric perturbations, whose invariance under fl56|) is again guaranteed by fll30p . 
(11311) . is given by fll28p . (11291) by virtue of (I126p . In this respect, these bound- 
ary conditions are the natural generalization of magnetic boundary conditions 
for Euclidean Maxwell theory, where one sets to zero at the boundary the 
tangential components of the potential, the gauge-averaging functional, and 
hence the ghost 0-form. The boundary conditions (I127p - (ll3ip were consid- 
ered in [36] as part of the effort to understand the relation between the wave 
function of the universe and the effective action in quantum field theory. The 
loop expansion in quantum cosmology was then obtained after a thorough 
study of boundary conditions for the propagator. 

In light of (I127p . the boundary conditions (I128p . (I129p lead to mixed 
boundary conditions on the metric perturbations which take the form (cf. 
[3811391110]) 



^^00 6 9 



0, 



(132) 



dM 



dhpi _^ 3^ _ 1 dhpo 
dr r °* 2 9x* 



0. 



(133) 



dM 



To evaluate the scaling behaviour of the corresponding 1-loop amplitudes, it 
is necessary to write down the mode-by-mode form of the boundary condi- 
tions ([I32D, ([I33D, ([HID, (II30D and ([131]). They lead to 



dan 6 
— h -ar, 

CLT T 



1 dCn 

dr 



dr 



-On 
T 

dCn 

dr 



(n^ - 1] 



at dM, 



at dM, 



+ 


3 

r 


= 


at dM, 


dn 


= 


at 


dM, 


Cn 


= 


at 


dM, 


fn 


= 


at 


dM, 


k 
rvn 


= 


at 


dM, 



(134) 

(135) 

(136) 

(137) 
(138) 
(139) 
(140) 
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/„ = at DM, (141) 

m„ = at dM, (142) 

p„ = at dM. (143) 

On using, for example, the technique of section 2, the corresponding contri- 
butions to C(0) are found to be (the results quoted below are independent of 
the particular algorithm used) 

278 

C(0)transverse—traceless modes . _ ; 

(144) 

45 

C(0) partially decoupled modes 2 15 17, (145) 

11 2 31 494 

C(0) vector modes — 12 — — ^ g ~ JgQ ~ ^5~' (146) 

15 

C(0) decoupled vector mode ^) 

(147) 

__ /179 ^\__149 

C (Uj scalar ghost modes I 120 360/ ,f cr ' I "J 



45 



_ / 41 31 \ _ 77 

C(0) vector ghost modes — ( ~ ~ 360 / ~ 90' (149) 

5 

C(0) decoupled ghost mode 7^' 

(150) 

7.1 Eigenvalue condition for scalar modes 

For scalar modes, which are not discussed in the previous list of results, one 
finds eventually the eigenvalues E = from the roots X of 



Ti 

j;(x) ± -Jnix) = 0, (151) 
X 

j;(x)+(-|±^) J„(x) = 0, (152) 

where J„ are the Bessel functions of first kind. Note that both x and —x 
solve the same equation. 
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7.2 Four spectral zeta-functions for scalar modes 

As we know from section 2, by virtue of the Cauchy theorem and of suitable 
rotations of integration contours in the complex plane [H |2T], the eigen- 
value conditions fllSip and (11521) give rise to the following four spectral zeta- 
functions 



^2s 



(153) 



where, denoting by J„ the modified Bessel functions of the first kind (here 

P+ =n,(3_ = n + 2), 

F^{zn) = z-^^ (^znl'^izn) ± nln{zn)^ , (154) 

F^{zn) = z-^± (^znl'^izn) + ± r?j /„(^n)^ . (155) 

Regularity at the origin is easily proved in the elliptic sectors, corresponding 
to Cl(s) and Cb(s) 



7.3 Regularity at the origin of ( 



B 



With the notation in Refs. [HI |T9], if one defines the variable r = (1 +z^) 2 , 
one can write the uniform asymptotic expansion of in the form 



On splitting the integral dr = dr + dr with yU small, one gets an 
asymptotic expansion of the left-hand side of Eq. f ll53p by writing, in the 
first interval on the right-hand side, 

°° r,- +(r) \ ^ Ri -i-(T) 



and then computing [HI [19] 

CAr) = ^ = (1 - rr-' E Kjf^r^. (158) 



OR, 



a=j-l 
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Remarkably, by virtue of the identity obeyed by the spectral coefficients 
on the 4-ball, i.e. 

4i 



a 



a=j 



T{a-j + l) ^ 



(159) 



which holds Vj = 1, oo, one finds [HI [19 



1 

hm sCUs) = ^ a(a - l)(a - 2)irf) = 0, 



(160) 



a=3 



and [T8l[T9] 



12 



, , 5 1079 Iv^ / x.^f^^ / 296 



4 240 2 



45 



(161) 



a=2 i=l 

where, on denoting here by ip the logarithmic derivative of the F-function 



ir(a + 1) 
QTia-2] 



- log(2) 



(6a2-9a + l)r(a-2) 



r(a + ll 



+ 2^(a + l) -^/'(a-2) -^(4) 
'-1V r . 2-j 



(162) 
(163) 



Equation fll59p achieves three goals: 

(i) Vanishing of the log(2) coefficient in fll6ip : 

(ii) Vanishing of YlT=i fU)9ij) in ([MD; 

(iii) Regularity at the origin of 

7.4 Interpretation of the result 

Since all other C(0) values for pure gravity obtained in the literature on 
the 4-ball are negative, the analysis here briefiy outlined shows that only 
fully diffeomorphism-invariant boundary conditions lead to a positive C(0) 
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value for pure gravity on the 4-ball, and hence only fully diffeomorphism- 
invariant boundary conditions lead to a vanishing cosmological wave function 
for vanishing 3-geometries at 1-loop level, at least on the Euclidean 4-ball. If 
the probabihstic interpretation is tenable for the whole universe, this means 
that the universe has vanishing probability of reaching the initial singularity 
at a = 0, which is therefore avoided by virtue of quantum effects [IHl ITU] , 
since the 1-loop wave function is proportional to a^^^-* [7|. 

8 The strong ellipticity issue 

The result outlined in section 7 is non-trivial because the zeta-function 
corresponds to the sector of the boundary-value problem for which strong 
elhpticity [26] fails to hold [391 SO]. We now define in detail this concept, and 
we are also going to discuss its relevance both for physics and mathematics. 

Let M be a smooth, compact Riemannian manifold endowed with a 
positive-definite metric g and assume that the boundary dM is smooth. Let 

V he a. vector bundle over M and C°°(y, M) be the space of smooth sections 
of the bundle V . With the introduction of a Hermitian metric E and the 
Riemannian volume element on M, the dual bundle V* can be identified with 

V and a natural inner product for the smooth sections of V can be defined. 
It is clear that the Hilbert space £^(V,M) is identified with the completion 
of C°°(V, M) with respect to the inner product. An operator of Laplace type 
is a map 

^ : C~(\/, M) C^{V, M) , (164) 

expressed as follows: 

^ = -g'^'^l^X + Q , (165) 

where denotes the connection on V and Q represents a self-adjoint en- 
domorphism of V . 

The boundary data for the Laplace operator under consideration can be 
written as 

^ ( ) . (166) 

where ip G C°°(V, M) and we have set 

^o(v5) = '^\dM , i'i{<^) = Vatv^Ism , (167) 
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with V N denoting the normal covariant derivative with respect to the bound- 
ary dM. By introducing the tangential differential operator on DM the 
boundary conditions can be written in a concise way as 

B^iljdv) = , (168) 

where the general form of By, which ensures the self-adjointness of the op- 
erator £, is [381 lag HO] 

«-=(A(i-n))' ^'''^ 

where 11 denotes a self-adjoint projector and A is a self-adjoint tangential 
differential operator satisfying the relation 

HA = An = , (170) 

and I is the identity endomorphism of V . The operator A can always be cast 
in the manifestly self-adjoint form 

A = (I - n) I ^ (rv. + v,r) + (i - n) , (in) 

where Vj represents the covariant tangential derivative, compatible with the 
induced metric on dM, and F* and S are endomorphisms satisfying the re- 
lations 

P = E-'^iVyE = -V\ S = E-^S^E = S , (172) 

nr = rn = = = . (173) 

Here, the bar denotes the adjoint in the space £^(V,M) and dagger is the 
Hermitian conjugate. 

It is instructive to notice that different choices for the projector 11 and the 
operator A lead to different types of boundary conditions. More precisely: for 
n = I and A = one obtains Dirichlet boundary conditions, for 11 = A = 
one recovers Neumann boundary conditions, and for 11 = and A 7^ 0, with 
A not a differential operator, the boundary conditions are reduced to the 
Robin type. 

The leading symbol of the Laplace operator 01651) is defined as 

ai(^;a;,0 = 1^1' ■ I = g''\x)Ub ■ I , (174) 
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where ^ G T*M is an arbitrary cotangent vector. For a non-singular Rie- 
mannian metric the leading symbol f ll7ip is positive-definite and, hence, the 
operator =Sf is elliptic [H]. In order to analyze the strong ellipticity condi- 
tion for the boundary-value problem fll65p and (11681) we need to introduce 
the leading symbol of the boundary- value operator B^. Such leading symbol, 
denoted by ag{B^), is defined as [10] 

'^»(«^)=(ff (i-n)) • (1^5) 

where, by exploiting the relation (I17ip . one has 

T = -^aL{A)=rp^, (176) 

with an arbitrary cotangent vector p on T*{dM), the cotangent bundle of 
the boundary of M. From the relations (11721) one can prove that the matrix 
T defined above is anti-self-adjoint, 

f = -T, (177) 

and it satisfies the equation 

UT = TU = . (178) 

In a neighbourhood of dM the Riemannian manifold M can be locally 
described by a direct product Q = [0, e] x dM. A local set of coordinates for fl 
is = x'^{r, x*) where r denotes the normal distance from the boundary and 
X* are the coordinates on the moved boundary dM{r) = {x G M \r{x) = r} 
with r G [0, e]. In this setting, the leading symbol of the Laplace operator ^ 
can be written as ctl^^; x,r, p,uj). Let us set r = 0, make the replacement 
u — 7- —idr, and consider the resulting differential equation 

[aL{^; X, 0, p, -idr) - A • I] 0(r) = , (179) 

with A G C — M"*", and whose solutions must satisfy the asymptotic condition 

lim 0(r) = . (180) 

r— >oo 

The boundary- value problem, consisting of the pair (11650 and (11680 . is 
said to be strongly elliptic with respect to the cone C — if, for any p G 
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T*{dM), A G C — R"*" and ip'^^, there exists a unique solution to the equation 
(11791) satisfying both the asymptotic condition (USUI) and the relation 

ag{Bj^){x,p)^y{<P)=^'^{<P) . (181) 

For an operator of Laplace type the differential equation (11791) can be 
written as 

[-d^ + \p\^ - X] 0(r) = . (182) 

The general solution to the above equation which satisfies the condition (I180p 
of decay as r — )■ oo is 

0(r) = xexp(-/ir) , (183) 

where we have set /x = ^/\p\'^ — A. The boundary data for the solution (I183P 
can be expressed as 



V'^(0) = 




(184) 



The boundary-value problem under consideration is strongly elliptic if the 
boundary data (I184p satisfy equation (I18ip . More precisely, strong ellipticity 
holds if the matrix associated with the linear system 

is invertible. One can show that the above system is equivalent to the set 

i^iI-^T)x = M-V^i- (186) 

Since the first equation in (I186p is independent of the second jlO], verifying 
the invertibility of the matrix in (I185p . and hence strong ellipticity, is equiv- 
alent to verifying the existence of a unique solution to the second equation 
in ( I186p for arbitrary ip^ and . A necessary and sufficient condition for the 
existence of a unique solution to the second equation in 0186p can be found 
to be 

det [pi -iT]y^O. (187) 

The matrix iT defined in (I176P is self-adjoint, and therefore its eigenvalues 
are real. It is clear that for any A G C — IR+ the quantity p = \/\p\'^ — A is 
complex and, hence, the matrix [pi — iT] is non- degenerate. For A G M_ the 
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quantity fi is real and satisfies the inequality fi > \p\. The above remarks 
together with the condition fll87p imply that 



\p\I-iT > . 



(188) 



By noticing that (|p|I — iT){\p\I + iT) = + we can conclude that the 
boundary-value problem (11651) and (I168P is strongly elliptic if the eigenvalues 
of the matrix are real and greater than — |pp, i.e. 



for any p G T*{dM). 

In the setting of Euclidean quantum gravity it has been shown [3S1 HO] 
that the eigenvalues of the matrix T are the following: 



Since the eigenvalues of the matrix are and — |pp the strong ellipticity 
condition (I189p in noi satisfied. This means, in particular, that the dynamical 
operator of the metric perturbations endowed with diffeomorphism-invariant 
boundary conditions in the de Bonder gauge is not strongly elliptic. 

9 Concluding remarks and open problems 

In order to avoid the problem of the lack of strong ellipticity in Euclidean 
quantum gravity one can consider various alternative approaches. 

When deriving the operator that describes the dynamics of the metric 
perturbations the particular choice of the de Bonder gauge renders the oper- 
ator of Laplace type but leads to the lack of strong ellipticity of the associated 
boundary-value problem. It remains to be seen whether strong ellipticity can 
be preserved if one considers instead dynamical operators on metric pertur- 
bations which are non-minimal. 

Another approach is to study Euclidean quantum gravity with non-local 
boundary conditions [12], or with boundary conditions which are not gauge 
invariant and that eliminate the occurrence of tangential derivatives [43]. 

Although viable, the alternative approaches mentioned above contain 
some problems and difficulties. It is, therefore, still unclear what is the 



Q{T^) = 0, 3?(I|pp + T2) > , 



(189) 




with degeneracy [m(m + l)/2 — 2] 
with degeneracy 1 
with degeneracy 1 . 



(190) 
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most appropriate way to solve the problem of the lack of strong ellipticity 
in Enclidcan qnantum gravity. The result of section 7, however, shows that 
there exists at least one background where a meaningful ({0) value is still 
obtainable despite the lack of strong ellipticity. At a deeper mathematical 
level, strong ellipticity makes it possible to define the heat operator, which 
is however not necessary in order to define the resolvent or complex powers 
of the given elliptic operator. For the latter two, one needs just a sector of 
the complex plane free of eigenvalues of the leading symbol. [We are grateful 
to Gerd Grubb for correspondence about this issue]. The investigation of 
other backgrounds might provide further examples of meaningful ({0) values 
despite violation of strong ellipticity, and their relevance for quantum cos- 
mology and/or quantum field theory should be assessed. 

Acknowledgments GE is indebted to Ivan Avramidi for scientific col- 
laboration on strong ellipticity in quantum field theory and quantum gravity, 
and is grateful to the Dipartimento di Scienze Fisiche of Federico II Univer- 
sity, Naples, for hospitality and support. KK would like to thank Stuart 
Dowker for the many years of fruitful and very enjoyable collaboration. Sev- 
eral of the papers written together and with others would never have seen the 
light of the day without his extremely valuable advice. The work of AK was 
partially supported by the RFBR grant No 11-02-00643. KK is supported 
by the National Science Foundation Grant PHY-0757791. 

References 

[1] Misner C W 1957 Rev. Mod. Phys. 29 497 

[2] Hartle J B and Hawking S W 1983 Phys. Rev. D 28 2960 

[3] Hawking S W 1984 Nud. Phys. B 239 257 

[4] Dowker J S and Critchley R 1976 Phys. Rev. D 13 3224 

[5] Hawking S W 1977 Commun. Math. Phys. 55 133 

[6] Henneaux M and Teitelboim C 1992 Quantization of Gauge Systems 
(Princeton: Princeton University Press) 

[7] Schleich K 1985 Phys. Rev. D 32 1889 



37 



[8] Barvinsky A O, Kamenshchik A Yu and Karmazin I P 1992 Ann. Phys. 
(N. Y.) 219 201 

[9] Barvinsky A O, Kamenshchik A Yu, Karmazin I P and Mishakov I V 
1992 Class. Quantum Grav. 9 L27 

[10] Kamenshchik A Yu and Mishakov I V 1992 Int J. Mod. Phys. A 7 3713 

[11] Moss I G and Poletti S 1990 Nud. Phys. B 341 155 

[12] Esposito G, Kamenshchik A Yu, Mishakov I V and PoUifrone G 1994 
Class. Quantum Grav. 11 2939 

[13] Esposito G, Kamenshchik A Yu, Mishakov I V and PoUifrone G 1995 
Phys. Rev. D 52 2183 

[14] Branson T P and Gilkey P B 1990 Commun. Part. Diff. Eq. 15 245 

[15] Vassilevich D V 1995 J. Math. Phys. 36 3174 

[16] Moss I G and Poletti S 1994 Phys. Lett. B 333 326 

[17] Esposito G, Kamenshchik A Yu, Mishakov I V and PoUifrone G 1995 
Phys. Rev. D 52 3457 

[18] Esposito G, Fucci G, Kamenshchik A Yu and Kirsten K 2005 Class. 
Quantum Grav. 22 957 

[19] Esposito G, Fucci, Kamenshchik A Yu and Kirsten K 2005 J. High 
Energy Phys. JHEP09 (2005) 063 

[20] Bordag M, Ehzalde E and Kirsten K 1996 J. Math. Phys. 37 895 

[21] Bordag M, Geyer B, Kirsten K and Ehzalde E 1996 Commun. Math. 
Phys. 179 215 

[22] Dowkcr J S 1996 Class. Quantum Grav. 13 585 

[23] Dowker J S 1996 Phys. Lett. B 366 89 

[24] Olver F W J 1974 Introduction to Asymptotics and Special Functions 
(New York and London: Academic) 



38 



[25] Thorne R C 1957 Philos. Trans. R. Soc. London 249 597 

[26] Kirsten K 2001 Spectral Functions in Mathematics and Physics (Boca 
Raton: CRC Press) 

[27] Elizalde E 1995 Ten Physical Applications of Spectral Zeta Functions, 
Lecture Notes in Physics m35 (Berlin: Springer- Vcrlag) 

[28] Bordag M, Klimchitskaya GL, Mohideen U and Mostepanenko VM 
2009 Advances in the Casimir Effect (Oxford: Oxford University Press) 
Kellogg O D 1954 Foundations of Potential Theory (New York: Dover) 

[29] Kellogg O D 1954 Foundations of Potential Theory (New York: Dover) 

[30] Kennedy G 1978 J. Phys. All L173 

[31] Moss 1 G 1989 Class. Quantum Grav. 6 759 

[32] Moss I G and Dowker J S 1989 Phys. Lett. B 229 261 

[33] Bordag M, Kirsten K and Dowker J S 1996 Commun. Math. Phys. 182 
371 

[34] Abramowitz M and Stegun I A 1970 Handbook of Mathematical Func- 
tions (New York: Dover) 

[35] Esposito G, Kamenshchik A Yu and PoUifrone G 1997 Euclidean Quan- 
tum Gravity on Manifolds with Boundary, Fundamental Theories of 
Physics 85 (Kluwer: Dordrecht) 

[36] Barvinsky A O 1987 Phys. Lett. B 195 344 

[37] DeWitt B S 2003 The Global Approach to Quantum Field Theory, In- 
ternational Series of Monographs on Physics 114 (Oxford: Clarendon 
Press) 

[38] McAvity D M and Osborn H 1991 Glass. Quantum Grav. 8 1445 
[39] Avramidi I G and Esposito G 1998 Glass. Quantum Grav. 15 1141 

[40] Avramidi 1 G and Esposito G 1999 Gommun. Math. Phys. 200 495 



39 



[41] Gilkey P B 1995 Invariance Theory, the Heat Equation and the Atiyah- 
Singer Index Theorem (Boca Raton: CRC Press) 

[42] Esposito G 1999 Class. Quantum Grav. 16 1113 

[43] Luckock H C 1991 J. Math. Phys. 32 1755 



40 



